Magnetic scaling in cuprate superconductors. 
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We determine the magnetic phase diagram for the YP^CuaOe+i and L^-iSr^CuCU systems 
from various NMR experiments. We discuss the possible interpretation of NMR and neutron scat- 
tering experiments in these systems in terms of both the non-linear cr-model of nearly localized spins 
and a nearly antiferromagnetic Fermi liquid description of magnetically coupled quasiparticles. We 
show for both the 2:1:4 and 1:2:3 systems that bulk properties, such as the spin susceptibiltiy, and 
probes at the antiferromagnetic wavevector (n, n), such as 63 Ti, the 63 Cm spin relaxation time, both 
display a crossover at a temperature T cr , which increases linearly with decreasing hole concentration, 
from a non-universal regime to a z — 1 scaling regime characterized by spin pseudogap behavior. We 
pursue the consequences of the ansatz that T cr corresponds to a fixed value of the antiferromagnetic 
correlation length, £, and show how this enables one to extract the magnitude and temperature 
dependence of £ from measurements of Ti alone. We show that like T cr , the temperature T* which 
marks a crossover at low temperatures from the 2 = 1 scaling regime to a quantum disordered 
regime, exhibits the same dependence on doping for the 2:1:4 and 1:2:3 systems, and so arrive at 
a unified description of magnetic behavior in the cuprates, in which the determining factor is the 
planar hole concentration. We apply our quantitative results for YBa2Cu307 to the recent neutron 
scattering experiments of Fong et al, and show that the spin excitation near 40meV measured by 
them corresponds to a spin gap excitation, which is overdamped in the normal state, but becomes 
visible in the superconducting state. 

PACS: 74.20.Mn, 75.40.Cx, 75.40.Gb, 76.60.-k 



I. INTRODUCTION 

The exotic magnetic properties of high-temperature superconductors have been extensively studied during the last 
few years, both because of their intrinsic interest and in the hope that these might provide insight in the physical 
origin of high-temperature superconductivity. Although considerable progress has been made in understanding the 
magnetic behavior of these materials, the basic issues of the doping dependence of the magnetic phase diagram and the 
relationship between transport properties and magnetic behavior in the normal state of the superconducting cuprates 
have not yet been completely resolved. In this paper we construct new magnetic phase diagrams for the YBa2Cu 3 06+ K 
and La2_ 2: Sr 2: Cu04 system on the basis of NMR experiments. Our approach makes evident the connection between 
the appearance of a spin pseudogap and magnetic scaling, and enables us to deduce from NMR experiments of the 
e3 Cu spin-lattice relaxation rate the magnitude and temperature dependence of the length, £, which characterizes the 
strength of the antiferromagnetic correlations found in both underdoped and fully doped cuprate superconductors. 

It is well established that—the undoped materials are described by the 2D Heisenberg model with a fairly large ex- 
change coupling J ~ 1550iul The thermal fluctuations destroy long-range order at finite temperature (the Hohenberg- 
Mermin- Wagner theorem) , and the Neel .transition observed in undoped cuprates is solely due to the small, but finite 
interplanar coupling. Chakravarty et ao gave a strong indication that the long- wavelength action of s = 1/2 2D 
Heisenberg antiferromagnet is reduced to the quantum non-linear sigma model, hereafter cr-model, characterized by 
an action, 



S, 



eff 



dr 



d 2 u 



(Vn) 



(0 T n) : 



(1) 



Here n is a three-component vector field, which describes the local staggered magnetization, which is subject to 
the constraint |n| = 1, p° s is the spin stiffness, and Co is the spin wave velocity. There is a short-distance cutoff 
A -1 for the spatial integrals. The properties of the cr-model are conveniently expressed in terms of the coupling 
constant g = Ticq/ p%. Chakravarty et aB identified three different regimes of behavior (Fig.|l|), defined by the zero- 
temperature ordered-non-ordered phase transition at g c = An/ A. According to the cr-model g — T phase diagram, for 
all coupling constants there exists a high-temperature Quantum Critical (QC) regime ((ps, A) < T < J) in which 
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the characteristic energy scale is set by temperature. At lower temperatures, depending on the coupling constant, g, 
the correlation length either grows exponentially, as in the Renormalized Classical (RC) regime (for g < g c ) or stays 
constant, as in the Quantum Disordered (QD) regime (for g > g c ). For the states without long-range order (g > g c 
or T > 0) the Goldstone mode is eliminated, which leads to a finite gap in the spin wave spectrum A = c/£, where c 
is the spin wave velocity, and £ the correlation length. 

Doping changes the ground state of the cuprates from a NeeirOrdered state to a highly-correlated ground state 
without long-range orderu. It was first argued by Sachdev and Yeu that the behavior of the cuprates for small doping 
is controlled by the zero-temperature critical point of the er-model leading to the dynamic critical exponent 2 = 1. This 
critical point does not necessarily correspond to the metal-insulator transition. The. Quantum Critical (QC) z = 1 
transition exhibited by the er-model was then analysed in detail by Chubukov et a£j, who emphasized its potential 
importance for the magnetic phase diagram of tho-,cupratc superconductors. n 

On the basis of NMR experiments of Imai et aEi and Takigawal, Sokol and Pinesu then argued that the cr-model, 
QC 2=1, description of the spin fluctuations continues to be valid at the comparatively large doping levels which 
are appropriate to the undcrdopcd cuprate superconductors. They suggested that a spin wave description is valid at 
high temperatures for these systems, while at higher doping levels, such as are found in YBa2Cu307, the particle-hole 
excitations act to completely damp the spin wave spectrum, so that the normal state of YBa2Cu307 is best described 
using the nearly antiferromagnetic Fermi liquid (hereafter NAFL) model which-iad been successfully applied to the 
analysis of NMR experiments in both its normal and superconducting statal3~ll3. In this mean field model the spin 
fluctuations at (%, 7r) are characterized by a relasational mode which varies as £ -2 ; thus it resembles a system with 
the dynamical critical exponent, 2 = 2. In Ref£3 the mean-field formula for the crossover from z — 1 to z — 2 was 
obtained by incorporating the quasiparticle damping in the cr-model expression for the spin susceptibility, and different 
experiments were analysed using this formula. A major step towards the description of the z — 1 to z — 2 crossover 
has been taken recently by Sachdev et ali-3, who explicitly derive the crossover scaling functions of the spin-fermion 
model using a 1/N expansion. 

A very, general analysis of both NAFL and nearly ferromagnetic Fermi liquid behavior has been carried out by Millis 
and IoffetSE3 while Monthoux and PineaLU, hereafter MP, have discussed the possibility that a NAFL approach, in 
which non-linear feedback from the magnetic interaction between quasiparticles acts to bring about a spin pseudogap 
in the quasiparticle spectrum, might yield the functional equivalent of QC z — 1 behavior over a broad temperature 
region in the underdoped cuprates, with overdamped spin waves arising because the magnetic energy, wj ~ a few 
J, and the quasiparticle Fermi energy are on approximately the same scale. They suggested that the onset of the 
pseudogap would occur at a certain value of the correlation length, which, is. independent of doping. 

In the present paper we extend previous analyses of NMR experimentsO'Ej on the metallic low temperature normal 
state regimes of La2- 3; Sr 3 ;Cu04 and YP^CusOe+z systems and present revised magnetic phase diagrams (Fig. || 
and Fig. |J) which make explicit the presence of both an upper limit, T cr , and a lower limit, T*, to QC 2 = 1 
scaling behavior in these systems. We show that T cr both marks the onset of-spin pseudogap behavior, a reduction 
of the effective quasiparticle density of states, found first in bulk susceptibility^ and Knight shiftily experiments, and 
subsequently in specific heat experimentsEll, and the crossover from a high-temperature non-universal regime to a 
low-temperature universal scaling regime. This scaling regime has dynamical critical exponent 2 = 1. We argue that 
this crossover happens at a fixed value of the correlation length, £ cr ~ 2, and that T cr is close to T x , the temperature 
which marks the maximum of the spin susceptibility. We use the results of NMR experiments to determine the region 
of applicability of scaling and find that T cr varies nearly linearly with doping level over a region which extends from the 
AF insulators, Yl^CusOe and La2Cu04, to YBa2Cu307 and the overdoped system Lai. 76 Sr .24CuO4. We find (see 
Fig. ^) that when the crossover temperature, T cr , is plotted as a function of planar hole doping, within experimental 
uncertainties, one finds the same dependence on doping for the 1:2:3 and 2:1:4 families. This demonstrates that the 
onset of magnetic scaling is determined only by hole concentration (band-structure effects are unimportant) and that 
bilayer coupling plays little or no role in determining spin pseudogap and scaling behavior. 

The QC-QD crossover at T*, which is predicted by the cr-model at lower temperatures, is also identified from the 
NMR data. Our identification of the experimentally determined "spin pseudogap" temperature, T cr , with a fixed value 
of £ then allows us to use NMR measurements of the 63 Cu spin-lattice relaxation rate to determine the correlation 
length and other parameters which characterize low frequency magnetic behavior as a function of temperature for 
both the YBa2Cu306+a; and La2- a; Sr 2 ;Cu04 families. We establish thereby the details of the variation with doping 
of the correlation length. 

There exists as well an intimate relationship between spin and charge response of the underdoped cuprate super- 
conductors. Ito et a£3 showed for YBa2Cu3 06.63 that below a temperature, T p , the planar resistivity ceases to exhibit 
a linear in T behavior, while the Hall effect likewise changes character at this temperature; a similar conclusion for 
YBa2Cu40g was reached by Bucher et acj. We find that within experimental error, T p = T*. In the La2- a; Sr 2 ;Cu04 
system Hwang et aM3 have shown that for x > 0.15, the Hall coefficient, Rh(T), follows scaling behavior, with a 
characteristic temperature, Th, which not only exhibits the same strong dependence on doping level that is found 
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for T cr , but which possesses very nearly the same magnitude as the T x determined from the maximum in the spin 
susceptibility. For this same system, with x < 0.14, Nakano et aO find that below T x , the resistivity ceases to be 
linear in T. Thus both the lower crossover temperature, T*, and the upper crossover temperature, T cr , which mark 
the limits of the QC z — 1 scaling behavior we identify in low frequency magnetic measurements, possess counterparts 
in transport measurements. 

Although there is indirect experimental evidence for their existence from NMR experiments, propagating spin waves 
for the underdoped cuprates have not been observed in neutron scattering experiments on the normal state. We show, 
on the basis of the parameters we infer from NMR experiments, that just such a result is to be expected - that 
spin waves are substantially overdamped by quasiparticle interactions in the normal state. Since this damping is 
markedly reduced in the superconducting state, we conclude that spin wave excitations should become visible in the 
superconducting state, and we ppesent arguments which suggest that this provides a quantitative explanation for the 
recent experiments of Fong et ac3 on inelastic neutron scattering from YBa2Cu307. 

Our paper is organized as follows. In Section II we review the one-component approach to the interpretation 
of NMR measurements. In Section III we use this to analyse the current experimental situation for the magnetic 
measurements in the underdoped YBa2Cu306+z and La2-xSr x Cu04 families, and on the basis of this analysis we 
arrive at our new magnetic phase diagram for these systems. In Section IV we discuss possible theoretical models, 
and the form for the spin susceptibility at high energy transfer. We discuss the application of the scaling analysis to 
inelastic neutron scattering experiments in section V, and in Section VI we present our conclusions. 



II. THE ONE-COMPONENT APPROACH TO NMR EXPERIMENTS. 



We start by giving a brief overview of the NMR experiments and their interpretation in one-component nearly 
antiferromagnetic Fermi liquid (NAFL) theory. The interested reader will find a moregdetailed discussion in the 
excellent recent review of NMR measurements in the cuprate superconductors by SlichterEI. 

Nuclear spins probe the local environment. The Knight shift provides a measure of the uniform magnetic suscepti- 
bility at a particular nuclear site, while the measurements of the spin-lattice relaxation rates yield information on the 
imaginary and real parts of the dynamic spin susceptibility x(q, The single-component description of thc-^JMR 
experiments in YBa2Cu307 has its basic justification in the observation by Alloul et aE3 and Takigawa et a£a that 
the 63 Cu, 17 O and 89 V Knight shifts see the same spin susceptibility. Since the 63 Cu spin-lattice relaxation time both, 
has anomalous temperature dependence and is much shorter than that for the 11 and 89 Y nuclei, it was proposedEZl 
that the single magnetic component formed by the system of planar_Cu 2+ spins and holes mainly resides on the planar 
copper sites. The Shastry-Mila-Rice (hereafter SMR) HamiltonianE3 which describes the coupling of the Cu 2+ spins 
to the various nuclei in YBa2Cu3C>7, has the form: 

nn 

Hmrs = 63 4(r 4 ) 5^A°%(r i )+B^5 a (r i ) + 

_ 3 

nn 

30 

nn 

+ S9 I a (r l )DY / S a (r J ), (2) 



where A a p is the tensor for the direct, on-site coupling of the 63 Cu nuclei to the Cu 2+ spins, B is the strength of 
the transferred hypcrfinc coupling of the e3 Cu nuclear spin to the four nearest neighbor copper spins, while O and 
89 Y nuclei see only their nearest neighbor Cu 2+ spins through the transferred hyperfine interactions C Q/3 and D. It 
follows from this Hamiltonian that different nuclei probe different regions in momentum space of x(q,,«|— > 0). Using 
eq. (H) , it is straightforward to express the spin contribution to the Knight shifts for the various nuclefcJ: 



= (A ll+ AB) X0 



6 3 K s _ ( A ± + 4£)XQ 



63,., -, * 2 



17 K S 



2^3X0 



~ 17 



3 



>K S = 



Here the j n are various nuclei gyromagnetic ratios, "f e is the electron gyromagnetic ratio, \o is the temperature- 
dependent static spin susceptibility, while the indices || and _L for copper nuclei refer to the direction of the applied 
static magnetic field with respect to the axis perpendicular to the Cu — O planes. The index (3 for the oxygen nuclei 
can be either || or _L The spin-lattice relaxation time, ( a Ti)p, for a nucleus a responding to a magnetic field in the 
j3 direction, is: 



a T- p \T) = 



2n 2 B h 2 uj ^ 



Q ^(q)x"(q,^^o), 



(4) 
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(5) 



The form factor 63 -F e // is the filter for the 63 Cu spin-echo decay time, 63 T 2 gS 



63-7^—2 



T 2 g(T) 



0.69 
1~28 



1/2 



l^F e// (q) X '(q) 



(6) 



The hyperfine coupling constants Ax, Ay, C, and D, which enter the SMR Hamiltonian, Eq. (^), were deter- 
mined using the.Kpight shift measurements in YBa2Cu30y and the frequency of the antiferromagnetic resonance 
in YBa2Cu306t3o. The momentum dependence of the form factors is shown in Fig. [| A striking difference of 
the tempepjbure dependence and the magnitude of the spin-lattice relaxation on the 63 Cu, 17 0, and 89 Y nuclei led 
Millis et al!3 to the conclusion that the magnetic response function is strongly peaked at the antiferromagnetic wave 
vector, as might be expected if one were close to an antiferromagnetic singularity. Since the form factors for 17 O 
and 89 Y nuclei vanish at Q = (it /a, n/a), the difference in the magnitude and temperature dependence of e3 Cu and 
other relaxation rates can be easily explained. Once the assumption of strong enhancement of the magnetic response 
function at the antiferromagnetic wave vector is made, the 63 Cu spin-lattice relaxation anisotropy ratio provides a 
self-consistency check on the values of the hyperfine constants, 



R = 



^(Q) 
3 W) 



The values of the hyperfine constants, 



B 



3.82 x lQ- 7 eV, 



A 



0.84B, Cu 



A,, = 
0.915, 



3.7 



-45, 



(7) 



(8) 



satisfy this check. Moreover, these hyperfine coupling constants remain unchanged (at the 57, to 107, level) with the 
variation of the doping level or for different cuprate materials, as might be expected from quantum chemical arguments. 
Thus these constants can be regarded as the same for both the YBa2Cu306+a: and I^-^Sr^CuC^ families. 

We follow Millis et aS3 and Monien et aSB and use a response function consisting of an anomalous part, Xmmp(<1, w), 
and a Fermi-liquid part, ^(q, oj), 



x(q, w ) = Xmmp (q, u) + xfl (q, w) . 



(9) 
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The expression used by Millis et a for the anomalous part can be written in the mean-field form: 

X(q, cj)mmp = r— 7 —7^ —, , (10) 

where 

XQ = < 2 (11) 

Here losf is the characteristic frequency of the spin fluctuations, £ is the correlation length, a is a scale factor. Due 
to the strong antiferromagnetic enhancement, the anomalous contribution to the copper relaxation rate dominates at 
all doping levels. However, since the form factors for other nuclei are zero at the antiferromagnetic wave vector, the 
Fermi liquid (Korringa-like) contribution to their spin-lattice relaxation rate becomes important. 

In the limit of long correlation lengths the expressions for the 63 Cu relaxation rates eqs (^), (|^) are considerably 
simplified^: 

63 (TiT) = 132 (a ■ K)/eV 2 u SF /a, (12) 

while 

63 J_ = 295 {eV/s) ai (13) 

J 2G 

Thus, up to a scale factor a, 63 (TiT) provides a direct measurement of u>sf, while 63 T 2 g provides a measurement of 
the correlation length. Moreover, the ratios: 

^J-=3.9x 10HK/eV)(to S F^ = c') (14) 

63^ 'J 1 

4- = 1.15 x 10 7 (K/s)(aujsF^ = Xq^sf) 



63^2 
± 1G 



(15) 



tell us-about the scaling lawsB, if any, obeyed by the low frequency magnetic excitations. Thus, as emphasized in 
RefsEl,E3, if 63 TiT/T 2 g is independent of temperature, one has QC z = 1 scaling behavior, and the doping-dependent 
ratio, ujsfS, — c', is determined. On the other hand, when 63 TiT/ 63 T| G = const, u>sf displays either non-universal 
mean field or QC z — 2 scaling behavior, and the product xq^sf is independent of temperature. 



III. EXPERIMENTAL PHASE DIAGRAMS. 



In this section we use the NMR results for 63 Ti, 63 T2G and the Knight shift on the underdoped 1-2-3 and the 2-1-4 
families to establish the complete experimental magnetic phase diagram for these systems. We first review briefly the 
properties of the regimes which might be relevant to the experimental situation in the metallic cupratesa: 

• The QC, z = 1 regime. In this regime the spin-spin correlator takes the form: 

X (q,^=£ 2 -"F(q£,^), (16) 

and the inverse correlation length varies linearly with temperature: l/£ = aT + b, where in the cr-modeli a and b 
are computable universal constants. The key signature of the QC z = 1 regime in NMR measurementscl is that the 
ratio 63 TiT / 63 T2G oc losfS, — d is temperature-independent (Eq. (|l~5|) ) , and c' is, for the cr-model, proportional to the 
spin-wave velocity c, d ~ 0.55c. In the QC regime the characteristic frequency losf is linear in temperature. The 
bulk magnetic susceptibility should also vary linearly with temperature, 

X(T) = Xumv(T) + xo, (17) 

where Xuniv(T) is a universal linear contribution, computable in the cr-modefl: 

Xuniv{T)= (^) 2 (0.34A :b T-0.137A(T = 0)), (18) 
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while the fermionic contribution, is temperature-independent. 



• The Quantum Critical, z = 2, regimeEjtJ. The characteristic energy scale in this regime is again set by temperature, 
so that at high enough temperatures ujsf oc T. However, since z = 2, ojsf oc £~ 2 . As a result, in this regime 
63 TiT/ 63 T| G = const, and at high temperatures the inverse squared correlation length should obey the Curie- Weiss 
law, £~ 2 = aT + b. 

• The QD regime. In this regime the correlation length (and 63 T2fj) saturates. If z — 2, the damping frequoaciz-ijs'f 
(and 63 TiT) should saturate at zero temperature to a certain value, determined by the fermionic dampingJlj'Eil. In 
case of the z = 1 cr-model, ujsf in the QD regime increases exponentially with decreasing temperatures. 

• Non-scaling regime. In this regime lattice effects are appreciable, and a scaling analysis is inapplicable. As a result, 
the calculations in this regime are model-dependent. A subclass of the non-scaling regime is the mean field NAFL at 
high doping, where Eq.(|To|) can be used. |— , 

In order to describe the z = 2 regime, the a- model Eq.([j]) should be generalized to include fermionic dampingta. 
The effect of this damping is that at low enough temperatures one always should find a crossover from the z = 1 
to z = 2 scaling regime. This crossover happens either in the QC or the QD regime, depending on the value of the 
fermionic damping. We briefly discuss this model in Section IV. 

We begin our analysis of magnetic behavior by considering La2-xSr x Cu04, where the spin susceptibility xo has been 
measured both in the bulk experimentsEfrtfLpid in NMR Knight shift probesE3cj, and 63 Ti measurements have been 
carried out for many different doping levelsQO. The decrease in the spin susceptibility as the temperature is lowered 
over a broad temperature region was the first identification in the experimental literature of the "spin pseudogap" 
phenomenon. At a temperature T x , the spin susceptibility (or the 63 Cu Knight shift) which at high temperatures 
typically decreases somewhat with increasing temperature, as shown in Figs and ^, changes behavior; below T x it 
decreases as temperature is lowered. Loram et a£B have carried out specific heat measurements which show that this 
decrease in xo is matched Jay a corresponding decrease in the quasiparticle density of states, N (T). Johnstonli3 and 
subsequently Nakano et aB3 also found that the bulk magnetic susceptibility could be expressed in a scaling form: 

X (T) = X o + n X (T/T max ) (19) 

Nakano et a@ used xoi K and T max as fitting parameters. They found that the temperature-independent term xo 
is doping-dependent; it therefore cannot be ascribed to the Van-Vleck contribution. This scaling behavior is easily 



explained if one identifies xo as the "fermionic" part, and K\(T '/T max ) as the "spin" part. Eq.(19) is then quite similar 
to Eq. (|l7|) , the result expected from the scaling theory. However, it should be emphasized that from our point of 
view, because of large lattice corrections, scaling ends above some temperature, T cr . To preserve Eq.(|l9|), the lattice 
corrections should enter "universally", i. e. preserving the relationship Eq.(|l7|), while breaking the cr-model linear 
temperature dependence of Xuniv(T)- 

Quite generally, for a z = 1 scaling regime which extends between a lower crossover temperature, T», and an upper 
crossover temperature, T cr , one expects scaling behavior to manifest itself at both long wavelengths (i. e. in Xo(T) 
) and at the wavevectors around Q which determine both 63 Ti and 63 T2G, with both Xo(T) and 63 TiT exhibiting a 
linear temperature dependence. T cr marks the onset of scaling behavior, while T* marks beginning of the crossover 
from the QC (z = 1) to thcj-QD regime, predicted in the scaling theoryEl, which has been previously identified in the 
NMR 63 Ti and 63 T 2G dataiM This is the temperature at which 63 TiT and 63 T 2G stop being linear in T, as is expected 
for the QC z = 1 regime. We propose that T cr marks the crossover from scaling and spin pseudo-gap behavior to 
non-universal behavior. Inspection of Fig.^Jthen show that T cr cannot be far from T x ; in other words, the maximum 
in Xo(T) is reached at very nearly the same temperature at which scaling begins. Inspection of Fig. |] shows that the 
crossover at T cr is also visible in 63 TiT, as a change in slope; for the various samples for which both Xo(T) and 63 7i 
measurements have been carried out, the two independent methods (one long wavelength, one short wave length) of 
determining T cr are seen to berin good agreement with one another. The lower crossover, at T», is readily visible in 
the 63 7i data of Ohsugi et ac3 and Imai et aE shown in Fig. [|. It is also clearly visible in Xo(T) measurements for 
Sr doping levels, 0.15 and below, as marking the onset of a more rapid fall-off in Xo(T) as the temperature is further 
decreased. Our experimentally determined values of T» and T cr are given in Fig. |jj|; T* is seen to vary comparatively 
slowly with hole concentration; it possesses a maximum for doping levels in the vicinity of Lai.ssSro.isCuO^ T cr 
on the other hand varies rapidly with hole concentration, for doping levels less than Lai.sSr .2CuO4. We indicate 
by a solid line our linear extrapolation of T cr to doping levels less than Lai.ssSro.isCuO^ (where high temperature 
measurements of 63 Ti have not yet been carried out); we note that this extrapolation suggests a crossover temperature 
of ~ 1200.fr for the AF insulator, La2Cu04, not far from that found from the high-temperature series studies of the 
2D Hcisenberg antiferromagnettj. 
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Spin pseudogap and scaling behavior has also been found in transport measurements on the-2-1-4 system. Haj 



ui- 



et aM5 and Nakano et 



2 A analysed the data for the bulk magnetic susceptibility, resistivitycj, and Hall effectE5E3 
measurements for this system. They identified T x from the spin susceptibility, and found similar characteristic 
temperatures, T p and Th, from an analysis of change with temperature of the resistivity and Hall effect. It is not 
necessary, from our point of view, that the temperatures identified using these criteria be exactly the same. However, 
as may be seen in Fig. ||, where our suggested doping dependence of T cr is compared with T x , T p , and Tr h , these 
respective temperatures are very close to each other; a closeness which makes evident the inseparability of spin and 
charge behavior in this system. We note that due to the vague nature of the crossovers large errors of ~ 157, are 
involved in our estimates for the values of crossover temperatures. Errors may also occur because the doping is not 
always exactly known, while different measurements are often performed on different samples. 

Having established that the phenomena of scaling and spin pseudogap behavior ace-, intimately related, and that 
both begin at T cr , we next explore the possibility suggested by Monthoux and PinestZI, that T cr is associated with 
a critical value of the correlation length, £ cr , which is independent of doping. One reaches the same conclusion on 
the basis of the cr-model, where one argues that lattice corrections bring about an end to scaling for £ >^ 1. We 
make the ansatz that £ cr ~ 2. Our reason for choosing this value is that one finds that in both the insulating and 
nearly fully-doped samples, the onset of spin pseudogap—behavior coincides with £ ~ 2. Thus in the insulator, the 
bulk susceptibility reaches its maximum at T ~ 1400idl3, not far from the temperature (T cr ~ 1000-K") at which 
£ = 2 according to high temperature series studiesEl, while, as we shall see below, at the opposite end of the doping 
scale, £ ~ 2 at the temperature (125 K) which marks the onset of spin pseudogap behavior in almost fully-doped 
YBa 2 Cu 3 07. 

Unfortunately, 63 T 2 g data is not yet available in the metallic 2:1:4 samples, so that one cannot verify directly (as 
we shall be able to do for their 1:2:3 counterparts) that the scaling law, lusfS, = c', a constant, is obeyed for these 
systems. However, on the assumption that it is valid , and moreover that £ cr = 2, we can determine the correlation 
length in the QC regime based on the NMR 63 T\ measurements only. To see this, note that between and T cr we 
have: 

63 7i(T)T = 132(s-K)-^—- (T*<T<T cr ), (20) 

so that a knowledge of T cr from another experiment (e. g. the maximum in xo{T)) serves to fix c'/a, and hence £(T). 
The doping dependence of c'/a can also be determined from 63 TiT only, since 

- = 63 (Ti(T cr )T cr )-^-- (21) 
a 66s • A 

The correlation length as a function of temperature and c'/a for different doping levels are shown in Figs |s|jic|. 

It is instructive to consider the extent to which the values of £(T) we have determined from 63 Ti display the 
universal behavior predicted by the cr-model for the z = 1 scaling regime, = a + bT, for the case in which c 
and c! do not vary with doping. We see from Table [j], the coefficients a and b are not universal, which is a clear 
indication that c, c', and hence a, must be doping-dependent in this system. On the other hand, we find that the 
slope of 63 TiT is nearly doping-independent between T* and T cr ( 63 TiT(ms ■ K) = 28 + 0.027T for Lai. 9 Sr .iCuO 4 
and 63 TiT(ms ■ K) = 10 + 0.025T for Lai.76Sro.24Cu04), which is what scaling theory for a doping-dependent c would 
predict. We note that our assumption that £ ~ 2 at T cr means that £ _1 may be written in a simple form, 



j =o+(i -a)(T/T cr ), (22) 

which experiments given in Table [l| satisfy. 

We turn next to the YBa2Cu306+a; family. A complete analysis can be done for the compounds YBa2Cu3 06.63j 
YBa2Cu408 and YBa2Cu307, where 63 Ti, 63 T2G, and Knight shift data are available. A less complete analysis, 
comparable to that given above for the 2-1-4 system, can be given for the Pr-doped YBa2Cu30y, where only 63 Ti and 
Knight shift experiments have been carried out. We start by analysing the Knight shift experiments in the Pr-doped 
compounds, Yi_ x Pr x Ba2Cu307E 3 J, with x=0.05, 0.1, and 0.15. We assume that Pr substitution acts to remove, on 
average, 1/2 hole per plane for each substituted Pr atom. Thus the planar hole concentration in Yi.g5Pro.o5Cu 3 07 is 
assumed to be 0.025 less than that of the "host" material, YBa2Cu307. But, as noted above, the hole concentration 
of that material is not precisely known, so that a corresponding degree of uncertainty affects the doping dependence 
we propose for this system. The spin part of the bulk magnetic susceptibility we obtain for these compounds (in 
states/eV) is shown in Fig.f?]. As we have done for the 2-1-4 system, we identify the onset of scaling behavior, T cr , 
with the maximum, T x , of the curve xo(T). We note that the bulk spin susceptibility varies linearly with temperature 
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between T» and T cr , as expected for the scaling regime. On turning to the 63 7\T plots, shown in Figs. |TT|a, [T^, we 
see that the lower crossover temperature, T*, is clearly visible, while, just as was the case for the 2-1-4 system, the 
crossover temperature, T cr , at which 63 TiT changes slope, is barely discernable. As we did for the 2-1-4 system, once 
T cr is identified, we can determine the doping dependence of both c'/a and l/£ for these systems, with the results 
given in Table [I] _ 

We next consider YBa2Cu307, for which measurements of the Knight shiftu, 63 Ti and 63 T2G have been carried 
out. The experimental results of Imai and SlichtciH for 63 T2G and 63 TiT in YBa2Cu3C>7 bring out two important 
points: 

• For T > 125 A ujsf and display Curie- Weiss (linear in T) behavior (Figs, [l^, |l4|) 

• For T > 125A cusfS, 2 = const. This behavior, which has been previously identified as the Gaussian, z — 2, scaling 
regime, could equally well be regarded as a manifestation of non-universal mean field behaviorc3. 

What happens at 125 A, where 63 TiT stops being linear in Tl The results of Walstedt et aS provide an essential clue. 
For their sample, as shown in Fig. [l3| the bulk spin susceptibility displays only a minor variation with temperature 
in YBa2Cu307, with a comparatively shallow maximum at 125K. However, at lower temperatures the variation of the 
spin susceptibility becomes much more pronounced. We therefore identify T cr as ~ 125A for the doping level present 
in this sample, which may be slightly on the underdoped side. Chir ansatz that £ = 2 at 125A then permits us to 
determine a directly from the 63 T2G measurements of Iff! et a£3; we obtain a = 15.6states/eV, a value slightly 
larger than that previously adapted by Thelen and Pineal, while, on making use of Eq. ^3, we find d = 35meV for 
this material. 

Our conclusion that T cr = 125A for "YBa2Cu307" thus assigns it a place as a full member of the 1:2:3 family, 
exhibiting the same magnetic scaling and spin pseudogap behavior as the underdoped system, Yo.95Pro.o5Ba2Cu30?, 
but with a substantially lower value of T cr . Viewed from this perspective, the 63 T2G measurements of Imai et a£3 
confirm our tentative hypothesis that the crossover at T cr is to a non-universal mean field regime. In order to examine 
more closely the doping dependence of T c ^Jpr this family, we need to know the oxygen doping level of the samples 
studied by Walstedt et a£3 and Imai et aSB. In the absence of detailed experimental information, we assume these 
correspond to YBa2Cu3 06.95, although a somewhat lower oxygen doping level would be equally plausible. We then 
obtain the results for the doping dependence of T cr and (c'/a) for the 1:2:3 family shown in Figs || and [l^. 

We next consider YBa2Cu40s, a material which possesses the same planar hole concentration as-YBa2Cu306.8j for 
which both NOR and NMR measurements of 63 Xi and 63 T2G have been carried out by Corey et acJ, with the results 
given in Fig. [Ll]a,b. They obtain the very important result that within experimental error the ratio, ( 63 XiT/T2g), 
is constant between 215A and 450A, thus demonstrating that YBa2Cu40s exhibits scaling behavior over a broad- 
temperature range, between T* = 215A and T cr > 450A. Our examination of the Knight shift results of ZimmermanEZl 
suggests to us that x(T) reaches its maximum value near 500.fr, so that T x ~ 500A; we therefore adopt the value, 
T cr = 470A", and use this assignment to obtain c'J-a — 4.57, from an extrapolation of the results for 63 Ti to this 
temperature. From the 63 T2G results of Corey et aEj, we find d = 50meV, and a = 10.9states/eV for this material, 
while the temperature variation of £ is given in Table | and Fig. 

The final member of the l:2r3 family for which 63 T2G measurements have been carried out is YBa2Cu30g.63j where 
the measurements of TakigawaO were used by Sokol and PinesH to establish scaling behavior, with T* ~ 170A. Because 
high temperature (T > 500A) measurements of 63 Ti 63 T2G or \ are not available for this system, one cannot obtain 
T cr from experiment. However, on assuming that T cr continues to vary linearly with oxygen doping, as Fig. ^| suggests 
that it does between YBa2Cu40s and YBa2Cu3 06.95, we find T cr = (840 ± 150) A, with the corresponding results for 
d , a, and l/£ given in Table | and Figs |l^, [l7]. Our results for the doping dependence of d and a show that while d 
decreases with increased oxygen content, a increases in such a way that the product ad is nearly constant, a result 
we examine further in the following section. 

Our determination of the dependence of T cr on oxygen doping for the 1:2:3 system makes clear the hierarchy of 
antiferromagnetic correlations shown in Fig. |l6|, in which ,at a given temperature, £ increases rapidly as the planar 
hole concentration is reduced. It also enables us to compare 1:2:3 system with the 2:1:4 system. To do so, we 
choose YBa2Cu306.5 as a basis, arguing that up to this oxygen concentration oxygen atoms added to the insulator, 
YBa2Cu306, go primarily into chain or interstitial positions, leaving the hole concentration in the plane near zero. We 
then assume one hole per plane for each added oxygen, so that YBa2Cu30g.95 corresponds to a hole concentration, 
rih = 0.225, etc, while for the La2- x Sr x Cu04 system we assume that each added Sr atom introduces a hole in 
the plane. The resulting dependence of T cr on hole concentration is given in Fig. ||. We see that within the not 
inconsiderable uncertainties of T cr and hole concentration, the two families display the same dependence of T cr on 
hole concentration. These results lead us to two important conclusions: 
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• Despite their somewhat different band structure and Fermi surfaces, the 2:1:4 system and 1:2:3 system possess 
essentially the same magnetic phase diagram. Put another way, that diagram is only very weakly dependent on band 
structure; it simply reflects the planar hole concentration. 

• The bilayer couplings!!] found in neutron scattering experiments on the 1:2:3 system plays little or no role in 
determining spin pseudogap and scaling behavior. 



IV. THEORETICAL MODEL 



Several models have been applied to describe the spin fluctuations in the metallic state. Two that appear consistent 
with experimental picture described above are the cr-modelaQ and the nearly antiferromagnetic Fermi liquid (NAFL) 
description The starting point of the er-model is a set of localized spins interacting with the fermionic background. It 
is possible then to integrate out fermions, so that one is left with the action of the ordering field. It can be shown 
that the a- model action Eq. (jl|) should then be generalized to include the effect of quasiparticle damping, if-this kind 
of damping is allowed by conservation laws. This leads to a description of the z — 1 to z — 2 crossoveiO. If the 
quasiparticle source of damping is absent, the fermions just change the cr-model coupling constant and the spin wave 
velocity. This model relies heavily on the existence of the propagating spin wave excitations in the metallic state, 
which so far has not been unambiguously confirmed in the direct experiments. In the NAFL description such spin 
wave excitations are not needed, but can arise as a result of a frequency-dependent effective interaction. The "spin 
pseudogap" is then related to a nonlinear feedback effect which brings about the equivalent of z = 1 scaling. We 
believe that both models describe from different viewpoints the same physics of the propagating spin excitations in the 
metallic underdoped cuprates. We now consider those in more detail, with particular attention to the generalization 
of the MMP model, Eq(Hj|), to finite frequencies, in order to be able to compare the model prediction with the results 
of neutron scattering experiments. 



A. QNLaM plus holes 

In this sub-section we discuss the spin-fermion model. We assume that the fermionic and spin degrees of freedom 
can be separated, and the action for the spins is the cr-model. We further assume a specific form for the interaction 
between the fermion and spin degrees of freedom: 

H int = Jo^T^+K-S)^, (23) 

where the ij are the spin indices. I— . 
It can be shownEa that integrating out fermions leads to the following action in this model-3: 

5 = ^ T E / 0|n(q,^)| 2 ( Co g 2 + ^+7Kl). (24) 

Here both g and cq are renormalizcd by fermions. Quasiparticles also lead to an additional damping, the "f\to n \ term 
in the action, if this kind of damping is allowed by the conservation laws. The er-model constraint \n\ = 1 remains 
valid in the this case. ._. 

Applicable to systems with two components^, the decomposition of the one-component system in two components, 
the "spins" S and the "fermions" ip could be artificial. In fact, at high doping YBaaCuaOg+z is known to be a 
one-component system, where the Fermi surface is measured in the ARPES experimentso. There is no proof that the 
modpL Eq. (|24|) applies in the one-component case, although some arguments in its favor have been given by Sachdev 
et AO. The necessary condition for the model Eq. (|24|) to be valid is , in our opinion, an almost local character of the 
copper spins. ._. 

The scaling analysis of the modeO specified by Eq.(|24|) leads to a magnetic susceptibility near the ordering wave 
vector (7T, 7r) which is a universal function: 

z ("he \* (q lu A T\ 

x(q '" )= (^r)^ \mT) *'[t>t>T>t)' (25) 



Here q marks the deviation from the ordering wave vector (it, ir). The two energy scales, A(T = 0) and T, reflect the. 
energy gap in the spin wave spectrum and the strength of the fermionic damping, respectively. The scaling analysistj 
shows that fermions do not necessarily overdamp and destroy the spin waves. At low damping the crossover is the 
same as in the pure cr-model. The only difference is in the presence of a finite fermionic damping in the QD regime. 
At higher T, there are two crossovers with decreasing temperature. The first one, from the high-temperature QC 
z = 1 to the low-temperature QC z — 2 regime, happens at T ~ T. At this temperature, according to this analysis, 
the temperature dependence of the correlation length £ cx 1/T should change to £ 2 oc 1/T. The lower second crossover 
at T ~ A 2 /r corresponds to the onset of the QD z — 2 regime. 

As is the case for the pure cr-modeltl, the a- model with fermionic damping, Eq. (|24|) , cap be approached using a 1 /N 
expansion technique. For N — oo, the resulting magnetic susceptibility takes the forrrJlj: 

at 2 

X(q, w) = — a (26) 

l + (q-Q)^ 2 -£ 

Here wgp ot £ -2 is the characteristic frequency of the damping in the particle-hole channel, A = c/£ is the value of 
the actual gap in the spectrum of the spin wave excitations of the disordered state^pEq. ( p(f ) can be understood as a 
the mean field propagator for the spin excitations damped in quasiparticle channelOcil. Note, that the only source of 
damping in mean field theory is the decay of spin wave into two fermions. Another source of damping, multiple spin 
wave scattering, appears only beyond the mean field, in the first order in l/A—Qn the contrary, the energy spectrum 
of the disordered magnons is already adequate in the mean field approximation^ 2 ]. To summarize these considerations, 
we therefore write: 

(27) 



^ F ermion ^Spin 

where uJF erm i on is the fermionic damping frequency appearing in the mean field theory, while tospin appears as a 
result of 1/N corrections, and displays spatial and frequency dependence on the scale of the correlation length. 



The imaginary part of eq.(26) has the form: 



x(q.">) = - J£ ¥ -> (28) 

from which we obtain the integrated (or local) magnetic susceptibility: 

hi \ f n 4. -ll^SF UJUJ SF 1\ 

Xl {uj) =na^- -tan \— — j J . (29) 



Note that for u> > A the imaginary part of the magnetic susceptibility Eq.(|2q) is peaked at an incommensurate 
wave vector, which is just a reflection of the spin waves at high energy. This fact is general and model- independent. 
Therefore, it should be seen either in a Fermi-liquid or a localized spin model. 

The spin wave excitations predicted by Eq.(|2^) need not be well-defined. A straightforward calculation shows that 
for q = Q, the maximum in x" occurs at 

uj max = A(l - ij), (30) 

where 

7 = A/lo sf , (31) 

from which one obtains: 

X"(q,aw) = ^ (32) 

7 

Hence for a well-defined spin wave (7 <C 1), one can find a resonant peak whose magnitude is far in excess of \Qi 
a situation which is to be contrasted with that which one obtains when only relaxational behavior is present, and 
Xmax = Xq/2- However, as we shall see, in the normal state of the cuprate superconductors, one likely has 7 > 1. 
Indeed, if we use the <T-model result, c' = 0.55c in the scaling regime, which extends down to < 100K in the 
systems studied in neutron scattering experiments, then, at (tt, 7t), 
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7=4- L8 ' T *^ T ^ T cr (33) 
c 

We note, that the relation that we have found above, ac' — const, follows directly from the sum rule for the localized 
spins, 

f ,2 [, X"{ci,u) , ,o„\ 

a q aoj 7=7 oc ac oc ac (34) 

J J l-e-"' 1 

As applied to the experimental situation described above, the spin-fermion model provides a reasonable qualitative 
description of the z = 1 scaling. Based on our analysis in Section III, we have not found signatures of the two 
crossovers in the scaling regime. From this, together with the results of the measured relaxation rate ratios, which 
yield z — 1, we would conclude that within the cx-model multiple spin wave scattering is the dominant contribution 
to damping of the spin waves in the underdoped materials. The detailed comparison of this theory with the inelastic 
neutron scattering experiments is then obscured by the fact that the correct form of damping for the spin wave 
excitations can only be obtained numerically for the experimentally accessible temperature and frequency regions. 
The " spin pseudogap" temperature in this model is the temperature where scaling stops working, i.e. the effects of 
the natural distance cutoff (~ the lattice constant, a) are appreciable. It is natural that this happens at certain value 
of correlation length, £ cr ~ 2a from Section III. Above the temperature T cr scaling analysis becomes inapplicable. It 
is impossible to obtain the precise value of T cr , or £ cr from the scaling analysis alone. A numerical estimate for £ cr in 
the insulator, as mentioned above, is consistent with ours. 



B. Nearly Antiferromagnetic Fermi Liquid Description. 

An alternative approach, the NAFL modf;]0, started from the high-doping level, where the magnetic fluctuations 
have been described as z = 2 GaussianEI The form eq. (|l0| ) is obtained by a Taylor expansion of the RPA 
susceptibility near the (n, 7r) point. The characteristic low-energy mode in this model is u>sf = r/£ 2 , with T of the 
order of Fermi energy. The spin waves are absent, and the spin fluctuations have rclaxational character. This model 
has been successfully applied to the explanation of the relaxation rates in YBa2Cu3C>7, as discussed in Section III. 

Recently, Monthoux and PinestH (hereafter MP) have proposed a phenomenological description of the observed 
spin pseudogap and z = 1 scaling behavior which begins at this "other end" , the high-doping limit in which a 
NAFL approach is the appropriate starting point. The NAFL description of magnetic scaling proposed by MP is 
based on the twin proposal that the strong Coulomb correlations between planar quasiparticles give rise to a second 
energy scale, uij ~ a few J, and that spin pseudogap behavior is a quasiparticle phenomenon, in which the dynamic 
irreducible particle-hole susceptibility, x(q, w), for wavevectors in the vicinity of (n, it) becomes temperature-dependent 
for T < T cr . More specifically, MP show that the propagating term in Eq.(p6|) can arise from the dependence on luj 
of the restoring force, / a (q, w, T), which gives rise to NAFL behavior in mean-field theory, in which case the spin gap 
is related to A by 

A = uj(xq/ X q) 1/2 =«>jilt (35) 

where xq = a t 2 is the nearly temperature-independent static irreducible particle-hole susceptibility. They also show 
that if one writes 

Lim^^oxiQ,^) = n q( t )xqu, (36) 

then in the mean-field approximation, 

USF = A^W' (3?) 
where £ ~ 7r _1 . Hence if the "commensurate" quasiparticle density of states, Nq(T), obeys the relations, 

Nq(T) = xq T > T cr (38) 

N Q( T ) = K 7 = ^77 T*<T<T cr , 
c 4 act, 

then one has the QC z = 1 scaling behavior, losf = c'/£, for T» < T < T cr , while above T cr , one recovers non-universal 
mean field behavior, with u>sf oc £~ 2 . 
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As MP show, since 



— ^- cx xqusf = Xq/Nq(T) (39) 

to the extent that Xq(T) is weakly temperature dependent, the temperature dependence of the quantity, Nq(T), may 
be directly obtained from experiment. Moreover, on making our ansatz, £ cr = 2, then at T cr we have, from Eq.(j39|), 
a scaling relation between a and c', 

ac' = I, (40) 

which is quite close to the values determined in Table In Fig |l8|, we give the results for Tf G / (T]T) for YIBa2Cu306.63, 
YB^C^Os, and YBa2Cu3C>7 obtained using Eq. (|39|) and the experimental results of Takigawa et au and Imai et 
aiEn. It is instructive to compare the temperature dependence of T 2 2 G /(TiT) oc (1/xq)Nq(T) between T* and T cr 
shown in Fig |lj for the YBa2Cu30g.63 an d YBa2Cu40s samples with that found for Xo(T), since at long wavelengths 
one expects to find, by analogy with Eq. (|36|) , 

Lim UJ ^ox"(Q,uj) = N (T)xo(T cr )(J 

-Xo(T)xo(T cr )^, (41) 

with N Q (T cr ) ~ xo(T C r) and N (T) cx Xo(T). We find that [^(TJ/xqCO] and Xo( T ) display a comparable temper- 
ature variation between T* and T cr ; however although both [Nq(T)/xq] and Xo(^) vary linearly with T over this 
temperature region, the respective slopes for the two quantities do not agree. 

We further note that the temperature which marks the lower endpoint of magnetic scaling behavior is barely 
visible in Fig. [l8| as the temperature at which [Nq(T)/xq] changes slope. From Eq.(|39|) we see that in the NAFL 
description of the QD state proposed in MP, at low temperatures, where £ ~ const, losf varies as (a + bT)^ 1 , as is 
found experimentally, rather than the exponential behavior predicted by the tr-model. 



V. INELASTIC NEUTRON SCATTERING RESULTS. 



In this section we discuss the consistency of the inelastic neutron scattering experiments with the results of NMR 
experiments and the scenario developed above. We start with YBa2Cu306+2;. Early experiments on this systemCj 
showed the presence of the gap in the excitations near the antifcrromagnetic wave vector, which increased with increas- 
ing doping, consistent with the scaling model discussed above. However, it has been difficult to extract the magnetic 
scattering from the unpolarized neutron experiments due to the presence of acoustic phonons at similar frequencies. 
The quality of the YBa2Cu306+ x samples and_measurements has improved considerably since the early work. In 
particular, an analysis of magnetic excitationscSEj has been done in YBa2Cu307 using a polarized neutron scattering 
technique, which enables one to extract the phonon part of the scattering cross-section. Detailed experimentso have 
also been performed on the underdoped compound YBa2Cu3C>6.6- In what follows, we focus on these two recent 
experiments.—. _ . 

Fong at a£El have recently shown experimentally that, contrary to the earlier work of Mook et no, the broad peak 
at 41meV corresponding to a magnetic excitation appears only in the superconducting state. As a result, they argue 
that this excitation indeed corresponds to quasiparticle pair creation in the superconducting state. Here we wish to 
suggest an alternative explanation for the 41meV magnetic peak: that it is a spin wave which only becomes visible as 
a distinct excitation in dhe superconducting state. Such an excitation is consistent with the singlet to triplet transition 
identified by Fong et ac3. As we have remarked in Section IV, the spin wave excitation is overdamped in the normal 
state. It is much better defined in the superconducting state, because NMR experiments show that lusf increases 
rapidly for T < T c ; indeed, for YBa2Cu30y, Martindale et a£3 find, ll>sf(0-75T c ) ~ losf{T c ) As a result, the gap 
in the spin excitation spectrum, A = e/£, should be observable in inelastic neutron scattering experiments. From our 
analysis for YBa2Cu30j we obtain d ~ 35meV and £(T C ) ~ 2.3. If A ~ AlmeV, then one must have 

c = A£ - 90meV (42) 

As a result, in the normal state, 7 = A/uj&tf ~ 2.5, and the spin wave excitations would be unobservable there, in 
agreement with the conclusion of Fong et a£3. On the other hand, at T ~ 0.75T C , one has 7 ~ 0.25, and the excitation 
becomes readily observable. 
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We next consider the cr-model and NAFL interpretation of this result. In the cr-model, if the damping of spin 
excitations is accounted by multiple spin scattering alone, the gap in the spectrum of magnetic excitations would be 
A(T C ) ~ 28meV. Since the actual gap is considerably larger, the fermionic source of damping must be appreciable 
in YBa2Cu3C>7. Using A = 1.82uspin£ an d Eq.(|27|), we obtain ^ s V i n I ^ f ermion — 0.47, as the ratio of the fermionic 
and spin sources of damping. Thus, the spin excitation in YBa2Cu30j could be a spin wave. As the temperature 
decreases below T c the particle- hole damping disappears. The gap in the spin wave spectrum causes the "spin" part 
of damping to also decrease with temperature. As a result, the spin wave spectrum is much better defined in the 
superconducting state. 

In the NAFL description, 

c = ujl (43) 

so that ujj ~ 2J, a not unreasonable value, while losf of course reflects only quasiparticle behavior, with its rapid 
increase below T c being entirely a consequence of the superconducting gap. Moreover, measurements of T^q in the 
superconducting state of YBa2Cu30st3 show that it gradually decreases below T c , a decrease which, according to 
Eq. (|l3|) , directly reflects a decrease in £. Since A = c/£, we would accordingly expect a corresponding increase in A. 
There are at present no measurements of Tig for YBa2Cu3C>7; if we assume that £ decreases, by some 20%, from 2.2 
at T c to 1.8 for T < T c /2 (by which temperature the gap is fully established), then a value of ~ 1.9 at 0.75T C is not 
unreasonable. On taking c = 74meV, we would then find A — 39meV at 70K, and A ~ AlmeV for T < 0.5T C , a 
result in reasonable agreement with experiment. 

We next estimate the spin wave spectrum in YBa 2 Cu3 06.55, where detailed inelastic neutron scattering results were 
reported recently by Sternlieb et aE3. From our analysis we find d ~ 70meV. In the cr-model, the "spin" source of 
damping may be expected to be dominant for this compound, and therefore c ~ 1277716^. Since this compound is close 
to the IM border, and exact doping level is not known, the correlation length is hard to determine by extrapolation. 
To obtain a gap of A ~ lOmeV, which is experimentally observed, we need £ ~ 12.7 at low temperature, a value 
which would require that the system be at a lower doping level than YBa2Cu30g.63- In the NAFL description, c 
would not be expected to change appreciably with doping, so that with c = lAmeV , we find a low temperature gap 
of lOmeV with £ = 7.4. This latter length is just what we find for YBa2Cu30g.63 near T c , so the NAFL description 
is clearly consistent with the emergence of a spin gap at lOmeV well below T c . 

Sternlieb et al found that the local (integrated) spin susceptibility exhibits ui/T scaling. The temperature 

at which the scaling breaks down decreases with increasing w. This scaling behavior is generic for the 2=1 QC 
regime. The scaling function for ^(cj/T) was found by Chubukov and SachdevQ. In the region of oj and T where the 
magnetic response was measured, x"(a->/T) oc lu/T, consistent with experiment. The decrease of the temperature T* 
with increasing energy transfer is also expected on the basis of scaling arguments. 

Provided 7 is small enough, the spin wave spectrum that we predict can be observed in neutron scattering exper- 
iments. According to Eq. (p6|), for uj > A x"(q, u>) is peaked at the incommensurate wave vector, displaced from 
(tt, 7t) by Sq = \J lo 2 — A 2 /c. However, direct measurement of the spin wave spectrum in YBa2Cu306.6 is obscured by 
some relatively large degree of imcommensuration observed in this materially. We estimate the width enhancement 
due to the spin waves at 30raeV for this material as 2Sq ~ 0.05[rJ.u.], while the experimental width of the peak 
is ~ 0.25 ± 0.1[r.l.u.] is independent of energy transfer. Thus, the predicted peak width enhancement agrees within 
experimental error with the value reported by Sternlieb et asB. ■— . 

We now consider Lai.ssSro.isCuC^, a compound where an extensive neutron scattering studyeJ of magnetic exci- 
tations has been reported. We note first that the magnetic excitations were found at incommensurate positions, in 
contradiction with the one-component theory of NMR discussed in Section II, which requires that x"(Qi w ) be peaked 
at a commensurate wave vector. According to Barzykin et so, the only possible resolution of this contradiction in the 
framework of the one-component model lies in the assumption of some sort of domain formation due to tiny regions of 
phase separation. If this is not what happens, more than one copper spin component should be involved. We further 
assume that the one-component theory is correct, and imply the presence of some internal superstructure leading to 
four distinct peaks in the neutron scattering experiment. However, we note that the analysis done here is more general 
than the one-component theory reviewed in Section II, since it is based on copper NMR only, and the possible incom- 
mensurability is not terribly important. Since 63 T2G experimental results are not so far available, to our knowledge, 
in the metallic Lai.ssSro.isCuO^ it is not possible to determine all the parameters entering theory We find that for 
a = 15States/eV, d ~ 29meV and £ ~ 7. Thus, this leads to a predicted spin gap, A ~ 7.5meV, for the cr-model, 
and A ~ lOmeV if c = 74meV for this compound. We also predict z = 1 scaling for T > 90K. This temperature T* 
should decrease with increasing energy transfer. None of these effects have been so far observed in Lai.ssSro.isCuO^ 
We found some vague indication of the spin gap of ~ 6meV in the integrated intensity measurements at 35K, where 
the most detailed experimental study was reported. However, further study at higher temperatures is required, in 
which the kind of scaling plot for the integrated intensity such as the one used by Sternlieb et al in YBa2Cu30g.6 
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would be extremely useful. We note that 2 = 1 scaling was found by Keimer et a in 2:1:4 system at lower doping 
level, in the insulating state, as expected from the model discussed above. 

The determination made here of the spin fluctuation parameters, a and £, for the 1-2-3 system enables us to estimate 
the maximum spin fluctuation signal strength ~ \q/2, measurable in inelastic neutron scattering experiments on the 
normal state. We find that \Q ranges from 69states/eV for YBa2Cri307 to 470states/eV for YBa2Cu306.63! this 
provides a partial explanation of the difficulty (compared to O6.63) in carrying out quantitative measurements on 
the fully-doped members of this family. As we have noted above, in the supercondicting state, once quasiparticle 
excitations can no longer damp the spin gap excitation at A, a considerably shorter signal, ~ {xq/i)i where 7 < 1/3, 
emerges from the continuum. We can carry out a complete estimate for xq f° r the 2:1:4 system if we assume that 
for this system as one hnds a result similar to that found for the 1:2:3 system, c'a ~ 0.54. We then find, through 
interpolation between our results for the SV0.13 and SV0.15 samples, that \Q ~ 790States/eV for Lai.ggSro.uCuC^ 
and xq ~ 2^00 States / eV for Lai.gSro.iCu04. thus one can anticipate significantly stronger spin fluctuation signals 
in the 2:1:4 system. It is interesting to note that the value of the correlation length at T* for La1.86Sro.14Cu.O4 is 
£ ~ 8, comparable to that found by Mason et asB in their neutron scattering experiments on this system. 

We do not attempt to compare, in this paper, the correlation length obtained in our analysis of the NMR exper- 
iments in YBa2Cu30g+a; system with that measured in the inelastic neutron scattering experiments. Part of the 
reason is that the observed peak at the commensurate wave vector is very broad, a result which would lead to an 
oxygen NMR spin lattice relaxation rate which is in disagreement with experiment. On the other hand, the neutron 
scattering experiments of Tranquada et ao in YBa2Cu30g.6 have suggested the presence of an unresolved four-peak 
superstructure. This superstructure, if present, would make such a comparison meaningless. 

VI. CONCLUSION. 

The "spin pseudogap" effect in the cuprate superconductors is a puzzling phenomenon which corresponds to a 
reduction of the effective density of states. We show that this effect is present, and has similar character, in both 
single-plane La 2 _ a; Sr a ,Cu04 and in bilayer YBa 2 Cu306+2; families. As a result, we have pursued a single-plane scenario 
for this effect in which any coupling between the bilayers does not play a crucial role. We find in this scenario that 
the development of the "spin pseudogap" corresponds to the onset of scaling at a certain value of correlation length. 
The value of the dynamical critical exponent in the scaling regime is z = 1. 

We have analysed a large collection of experiments using the scaling theory (the QNLcr-model) and the NAFL 
approach The scaling theory is based on the existence of spin wave excitations in the metallic cuprates, which have 
a finite energy gap at (71% w) and a linear dispersion law at high frequency. It should be emphasized, that there is 
no gap other than the gap of the spin wave excitations in the er-model. This gap should be observable in inelastic 
neutron scattering experiments in the superconducting state, and may have been seen as the 41 meV excitation in 
YBa2Cu30y. The "spin pseudogap" observed in the bulk magnetic susceptibility or specific heat measurements is of 
different character, it does not necessarily correspond to a gap in the quasiparticle spectrum. We relate the observed 
spin pseudogap behavior to the onset of scaling behavior at high temperatures. It should be noted that a scenario in 
which the spin pseudogap phenomenon is related to a failed SDW gapEil is also possible. This mechanism is somewhat 
similar to the NAFL framework, considered in Section IV, in which novel physics arises if one takes into account 
vertex corrections brought about by stroing magnetic correlations. 

In summary, we have applied the z = 1 Nearly Antifcrromagnctic Fermi Liquid Scaling theory to the analysis 
of the bulk susceptibility, NMR 63 Ti and 63 T2G relaxation rates, specific heat, resistivity, Hall effect, and inelastic 
neutron scattering experiments in the La2- ;r Sr :r Cu04 and YBa2Cu30g+a: families. The magnetic phase diagram for 
the two compounds display the same behavior, and can be superposed on one another, when plotted as a function 
of hole concentration. Two magnetic crossovers are evident from experimental data. The upper crossover, at T cr , 
corresponds to the onset of the z = 1 scaling behavior. We argue that above T cr the spin dynamics is non-universal, 
i.e. influenced by the lattice cutoff. As a result, no definite predictions can be made on the basis of scaling theory for 
T > T cr . Below T cr scaling applies, and universal scaling functions can be computed. The temperature T cr decreases 
with increasing doping. The lower crossover temperature , T* , marks the low-temperature end of the z = 1 Quantum 
Critical scaling behavior, and its appearance is predicted by the scaling theory. The magnetic susceptibility, the NMR 
relaxation rates 63 TiT, 63 X2Gj an d the resistivity are linear in the Quantum Critical regime between T* and T cr . Both 
crossovers disappear in the overdoped case, as the correlation length becomes too short for z = 1 scaling to apply. 
These considerations are summarized on Figs. |^, ||. There is no region of applicability for z = 2 Fermi Liquid scaling 
on our magnetic phase diagram in the normal state. The observed ratio, 63 TiT '/ T%q = const in YBa2Cu307, in our 
interpretation, corresponds to a mean field non-universal behavior. 

On making the ansatz that £ cr = 2 at the critical temperature T cr which marks the onset of the z = 1 QC 
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regime, we have determined, from the NMR experiments, the correlation length as a function of temperature for both 
YBa2Cu306+£i; and La 2 _ ;c Sr ;c Cu04 families, and the spin wave velocity in YBa2Cu306+a; as a function of doping. 
From this analysis, we can approximately obtain the energy spectrum of the spin waves, which should be observed in 
inelastic neutron scattering experiments on the superconducting state. Thus, in our scenario, the peak in magnetic 
inelastic neutron scattering in YBa 2 Cu307 at AlmeV, and the spin gap of WmeV observed in YBa 2 Cu 3 06.6 could 
originate in spin waves. The observation of the full spin wave spectrum is obscured in these materials by strong 
phonon scattering and by some degree of discommensuration found in YBa2Cu3C>6.6- 

In arriving at our magnetic phase diagram, we have frequently had to resort to interpolation and extrapolation 
from existing experimental data. Our calculations are, in many cases, directly subject to experimental test. It is 
our hope and expectation that the NMR experimental community will now fill in the interstices, and so make more 
precise the dependence of T cr and T* on hole concentration. For example, detailed measurements of 63 T cr , 63 7\ and 
Xo(T) carried out on the same samples, for samples of known oxygen concentration in the vicinity of YBa2Cu3 06.95, 
can not only determine T cr in this region, and establish its dependence on hole concentration, but can also verify 
whether, as we have conjcctur, the ratio, 63 TiT/ 63 T2G, becomes independent of temperature below T cr . In similar 
vein more detailed measurements on the Knight shift and 63 Ti of the overdoped 2:1:4 system will help one to deduce 
T* and T cr for this system. 

We are especially grateful to Takashi Imai, Charlie Slichter, and Bob Corey for communicating their data prior to 
publication and for numerous discussions of the above topics. We also thank Lev Gor'kov, A. Millis, Q. Si, A. Sokol, 
B. Stojkovich and Y. Zha for the discussions of these and related subjects. This work was supported by NSF grant 
No DMR9 1-20000 through the Science and Technology Center for Superconductivity. 
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FIG. 1. The phase diagram of the cr-model - from Chakravarty et aE 

FIG. 2. The magnetic phase diagram determined from NMR experiments on the YBa2Cu306+j/ system. The temperature 
T cr marks the crossover from non-universal behavior to spin pseudogap behavior, while T* marks the crossover from QC(z = 1) 
scaling behavior to QD behavior. Scaling behavior is thus found between T* and T cr . 

FIG. 3. The magnetic phase diagram determined from NMR and bulk susceptibility experiments on the La2- :E Sr a ;Cu04 
system. The solid circles represent values of T cr determined rby us from NMR experiments; the triangles correspond to the 
maximum-in the bulk susceptibility reported by Hwang et a£3; the open squares are the corresponding results obtained by 
JohnstorJia 

FIG. 4. The dependence on hole concentration of the temperature, T cr , for the onset of spin pseudogap behavior determined 
here for the 2:1:4 and 1:2:3 systems is compared to the characteristic temperatures determined by Hwang et aEB from transport 
and bulk susceptibility measurements and the maximum in bulk susceptibility found by Johnstonc3. 

FIG. 5. Form factors as a function of momentum, for planar oxygen and copper sites, in units of the hyperfine coupling 
constant, B 2 

FIG. 6. T cr (denoted by arrows) for Lai.gSro.2Cu04 and Lai.76Sro.24Cu04, es determined from the endpoint of linear in T 
behavior and the maximum of the 63 Cu Knight shift measured by Ohsugi et a£3J. 

FIG. 7. The spin susceptibility inferred from the Knight shift data for YBa2Cu3 06.63, YBa2Cu40g, YBa2Cu307, and 
Yi-zPr^BaaCuaOy, with x=0.05,0.1,0.15. 

FIG. 8. The temperature T* (denoted by arrows) for the QC-QD crossover, as determined from the 63 T\T measurements of 
Ohsugi et a£j (a) and Imai et all (b) in La2-xSra;Cu04. 
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FIG. 9. The temperature dependence of the correlation length in the QC regime, obtained by us for the L^-^Sr^CuO^ 
system. 

FIG. 10. The variation of c'/a with doping in La2_ I Sr a; Cu04 

FIG. 11. The determination of T„, the QC-QD crossover, for the three compounds YBa2Cu306.63, YBa2Cu408 and 
YBa2Cu30?: (a) From 63 Ti measurements (b) From 63 T2G measurements 

FIG. 12. The determination of T„, the QC-QD crossover, from 63 7i measurements on Yi_ a; Pr 2 ,Ba2Cu307. 

FIG. 13. The crossover behavior at 125K displayed by the spin contribution to the magnetic susceptibility!! and 63 TiTEl m 
YBa 2 Cu 3 07. 

FIG. 14. The Curie- Weiss temperature dependence for £ 2 determined from the 63 T2G measurements of Imai et aE in 
YBa 2 Cu 3 07. 

FIG. 15. c'/a, determined from the 63 Ti measurements, as a function of doping in YBa2Cu307_a. 

FIG. 16. The antiferromagnetic correlation length in the YBa2Cu307_i system, as determined from the 63 T2G measurements. 
For the Pr-doped compound (the line) the correlation length in the QC regime was determined from the 63 J\T measurements. 

FIG. 17. The parameters a (a) and c' (b), determined for the three compounds YBa2Cu30s.63, YBa2Cu40s, and 
YBa2Cu407, where both 83 Ti and 83 T2G relaxation time measurements are available. 

FIG. 18. 63 (T2 2 G /riT) (oc Nq(T)/xq), extracted from NMR experiments on YBa 2 Cu3 6 .63, YBa 2 Cu 4 08, and YBa 2 Cu 3 07 
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TABLE I. Spin fluctuation parameters determined from fits to NMR experiments. 



System 


T cr (K) 


T*{K) 


c'/a 
(10" 3 eV 2 ) 


Q 

(states /eV) 


c' 

(meV) 


ac' 


1/C (T. < T < T cr ) 


YBa9CmO(; <™ 


840 


170 


7.67 


8.34 


64 


0.53 


0.052 + 5.93 x 10~ 4 T 


YBa2Cu4 0g 


470 


215 


4.57 


10.9 


50 


0.545 


0.068 + 9.2 x 10~ 4 T 


YBaoCu-^O? 


125 




2.23 


15.6 


35 


0.545 




Yn Q^Pfn n^RfloChl*-* O7 
j-u.yo- 1 i u.uo i-j<^z Li o v j ( 


200 


115 


2.83 








0.32 + 9.0 x 10~ 4 T 


Yo.9Pro.iBa2Cu307 


2f5 


125 


3.0 








0.31 + 8.8 x 10~ 4 T 


Yo.85Pro.i5Ba 2 Cu 3 07 


260 


130 


3.18 








0.29 + 8.1 x 10~ 4 T 


Lai.9Sro.iCu04 


680 


50 


2.7 








0.024 + 7.01 x 10~ 4 T 


Lai. 85 Sro.i5Cu04 


4f0 


70 


1.91 








0.065 + 1.06 x 10~ 3 T 


Lai.8Sro.2Cu04 


f20 


65 


0.86 








0.23 + 2.25 x 10 _3 T 


Lai. 76 Sro.24CuC>4 


too 


35 


0.83 








0.255 + 2.46 x 10~ 3 T 
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